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Abstract 

We discuss the existence and characterization of quasi-stationary distributions and Yaglom 
limits of self-similar Markov processes that reach in finite time. By Yaglom limit, we mean 
the existence of a deterministic function g and a non-trivial probability measure v such that 
the process rescaled by g and conditioned on non-extinction converges in distribution towards 
v. If the study of quasi-stationary distributions is easy and follows mainly from a previous 
result by Bertoin and Yor [5] and Berg [2] , that of Yaglom limits is more challenging. We will 
see that a Yaglom limit exits if and only if the extinction time at of the process is in the 
domain of attraction of an extreme law and we will then treat separately three cases, according 
whether the extinction time is in the domain of attraction of a Gumbel law, a Weibull law or 
a Frechet law. In each of these cases, necessary and sufficient conditions on the parameters 
of the underlying Levy process are given for the extinction time to be in the required domain 
of attraction. The limit of the process conditioned to be positive is then characterized by a 
multiplicative equation which is connected to a factorization of the exponential distribution in 
the Gumbel case, a factorization of a Beta distribution in the Weibull case and a factorization 
of a Pareto distribution in the Frechet case. 

This approach relics partly on results on the tail distribution of the extinction time, which is 
known to be distributed as the exponential integral of a Levy process. In that aim, new results 
on such tail distributions are given, which may be of independent interest. 

Last, we present applications of the Frechet case to a family of Ornstein-Uhlenbeck processes. 

AMS subject classifications: 60G18; 60F05; 60G51. 
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1 Introduction and main results 

Let X be a continuous time strong Markov process with values in M + . We will denote by its 
distribution started at x and F u its distribution when Xq is distributed according to a probability 
measure v. It is supposed that there exists some a > such that this process is self-similar with 
index 1/a > 0, that is for c > 0, 

the distribution of {cX tc - a , t > 0} under F x is F cx for all x ^ 0. 

This defines a 1/a-positive self-similar Markov process (1/a-pssMp). We will assume furthermore 
that 

T := inf{u > : X u = 0} < oo F x -a.s. for all x > 0, 
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and that is a cemetery state for X. We denote by tp the supremum of the support of this extinction 
time under Pi, that is 

t F :=sup{i > : Pi(T < t) < 1}. (1) 
The goal of the paper is then twofold: 

• investigate the quasi-stationary (QS) distributions of X, that is study the existence and 
characterize the probability measures \i on (0, oo) such that for all t £ [0,tp), 

F^(X t GB\t< T ) = P M (B) 

for all Borel sets B; 

• study the asymptotic behavior of Xt conditioned on non-extinction. More precisely, our goal 
is to determine under which conditions there exists a function g : [0, oo) — > (0, oo) and a 
non-degenerate probability measure v on (0, oo) such that 

in) ( -^t | rp \ weakly 1 
Pi —TT £ * * < -M) > V 

\g(t) 1 °J t^t F 

We will then say that v is a limit in the Yaglom sense of X. It is not hard to see that the 
function g and the measure v are then uniquely determined up to a scaling constant, in the 
sense that if a convergence of type ([2]) holds for two pairs (g, v) and (</ , z/), then g(t)/g'(t) — > c 
as t — > ip for some c > and v' = v o (cid) -1 . 

In related literature, the terminology Yaglom limit is often used for convergences of the type ([2]) 
with g(t) = 1, t > 0. As we will see, for self-similar Markov processes, the limiting distribution of 
Xt conditioned on non-extinction is often trivial, converging either to the Dirac measure at or to 
the Dirac measure at 00. To evaluate the asymptotic behavior of Xt conditioned on non-extinction 
more accurately, we therefore have to normalize the process appropriately. 

At this stage we would like to mention that when ip = 00, if one replaces the definition ([2]) of 
Yaglom limit by asking that there exists a function g : [0, 00) —> (0, 00) such that for all x > 

P* G • I t < T ) ► 

\g(t) J t^oc 

with v x a non-degenerate probability measure on (0, 00) then, by the self-similarity of X, the 
random variable xXt/g(tx a ) conditioned on non-extinction has a non-trivial limit in distribution, 
and therefore, g(tx a )/(xg(t)) — > c x for some c x > 0, Vx > 0. This implies that g is regularly varying 
at 00. But, as we will see further, convergence in the sense ([2]) may hold for some functions g that 
are not regularly varying. Hence the definition ([3]) is more restrictive than ([2|). 

The problem of existence of QS distributions and Yaglom limits of certain classes of pssMp has 
been considered by some authors. For instance, Kyprianou and Pardo studied the case of a 
stable continuous state branching processes (see also Lambert [27] for similar results on non-stable 
continuous state branching processes). In [25], the first author considered the case of non-increasing 
pssMp (with some additional hypotheses) in order to describe the asymptotic behavior of solutions 
to some fragmentation equations. Our main purpose here is to provide a unified and general 
approach for QS distributions and Yaglom limits of pssMp that die at their first hitting time at 0. 
Specific applications will then be discussed, for instance to stable processes killed at (—00, 0], Bessel 
processes killed at their first passage at 0, stable continuous state branching processes, etc. Our 
main results can also be used to study the existence of quasi-stationary distributions and Yaglom 

1 In all the paper, convergences of probability measures are weak convergences. This is implicit from now on and 
will not be mentioned further. 
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(2) 



(3) 



limits for processes of the Ornstein-Uhlenbeck type associated to pssMp, particular cases of which 
are here studied. In particular in Example [3] we obtain a result closely related to Mandl's [29] 
seminal result about Yaglom limits for the classical Ornstein-Uhlenbeck process . 

In order to state our main results we recall first a few facts about self-similar Markov processes, 
Levy processes and exponential functionals of Levy processes. Further details will be given in 
Section 12.11 

It is well-known (see Lamperti [28] ) that there exists a KU{- oo} valued Levy process £ (with 
— oo as a cemetery state) independent of the starting point Xq, such that 

X t = X e W ^ T{txr) ), t>0, (4) 



where r is the time-change 

r(t) 



inf|s>0: J exp(a£ n )du > i j , t > 0, 



with the usual convention inf{0} = oo. The lifetime or first hitting time of — oo for the Levy process 
£ associated to X via Lamperti's transformation will be denoted by £. Lamperti proved that, for 
any x > 0, Tq is finite P x -a.s. if and only if either £ < oo a.s., or £ = oo a.s. and lim^oo £ 4 = — oo 
a.s., which are the cases we will generally consider in this paper. Conversely, Lamperti proved that 
given a Levy process £, the transformation just described gives rise to a 1/a-pssMp. We will refer 
to this transformation as Lamperti's transformation. Throughout this paper P will be the reference 
measure, and under P, X will be a pssMp and £ the Levy process associated to it via Lamperti's 
transformation; as already mentioned, P x then denotes the law of X issued from Xq = x. 

It follows from Lamperti's transformation that under P x the first hitting time of for X, Tq 
has the same law as x a exp(a£ s )ds. The random variable / defined by 

I := [ exp(a£ s )ds, (5) 



J o 

is usually named exponential functional of the Levy process £. Lamperti's above mentioned result 
implies that / is a.s. finite if and only if £ < oo a.s. or £ = oo and lim^oo^ = — oo, a.s. Using 
this it is easy to see that tp, as defined in (TjQ), is finite if and only if — £ is a subordinator with a 
strictly positive drift and that necessarily the drift is then equal to l/(atp) (see Lemma l2.ip . 

In the theory of self-similar Markov processes it is well known that exponential functionals of 
Levy processes are keystone, for instance they are used to describe entrance laws and allow to 
determine the local and asymptotic behavior of pssMp, see for instance [1] , [6] , [15] , [37] , [38] . So, it 
should not be surprising that they also play a central role in our development, nevertheless it is 
rather unexpected that the study of the maximal domain of attraction of exponential functionals 
of Levy processes allows to determine the existence of Yaglom limits, as it will be seen below. In 
section |2~TT1 we will provide further details about exponential functionals of Levy processes. 

For notational convenience we assume from now on that the self-similarity index is 
a — 1 (when X is a 1/a-pssMp, the process X a is a 1-pssMp and hence the transitions from 
1-pssMp to 1/a-pssMp are straightforward: when £ is the Levy process related to X via Lamperti's 
transformation Q, q£ is the one related to X a ). 

We recall that in the case where — £ is a subordinator, that is it has no-decreasing paths, its 
distribution is characterized by its Laplace exponent, </>, which is given by 

-ln(E[exp(A£ t )]) = t</>{\) 

(6) 



t(q + d\ + (1 - e- Xx )U(dx)^j , A>0, 
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where q > 0, d > and / °°(1 A x)U(dx) < 



oo. 



It turns out that the existence and description of QS distributions of pssMp is easy to settle. 
It mainly relies on results by Bertoin and Yor [5] who proved that when — £ is a possibly killed 
subordinator, there exists a random variable R independent of I such that 

RI L = e, (7) 

where e has an exponential distribution with parameter 1. The distribution of R is then uniquely 
determined by this equation and is a QS distribution of X. Its entire moments can be expressed 
in terms of the Laplace exponent <j> of — £ : 

n 

E[R n ] = Y\ VraGN. (8) 

1=1 

We denote here by fj,j the distribution of this random variable R. The results of Bertoin and Yor 
can be completed as follows. 

Theorem 1.1. X admits a QS law if and only if X has non-increasing paths. More precisely, 

(i) if X has non-increasing paths, the set of QS distributions of X is exactly 

{ / u5 e) o(Aid)- 1 ,A>0}; 

(ii) conversely, if X admits QS distributions, there exists a random variable R independent of I 
such that RI L = w e, which then implies that — £ is a subordinator. 



The proof of this result will be given in Section 12.21 

Studying the large times behavior of Xt conditioned on non-extinction is more challenging. A 
simple but key result for our analysis is the following equality of measures: for all t > 0, 

P (J - t E • | I > t) = Pi (xj £ • | t < T ) , (9) 



where / has the same law as / and is independent of (X s ,0 < s <t). This identity is proved in 
Lemma 12.21 Existence of a Yaglom limit for X is therefore connected to the behavior as t — > tp of 
the residual lifetime at t of the random variable /, more precisely it is equivalent to the existence 
of a function g such that the weak limit 

lim P ( — -— G • \ I > t ) exists and is non-trivial. 



According to the Pickands-Balkema-de Haan Theorem [TJ this, in turn, holds if and only if 
/ is in the domain of attraction of an extreme value distribution. The three families of extreme 
value distributions being the Gumbel family (we then write that I £ MDAQ um b e i), Weibull family 
(J G MDAwcibuii) and Frechet family (/ G MDAp^chet)- Some background and references on 
extreme value theory are given in Section 12.31 where the above assertions are also settled properly. 
See in particular Proposition 12.71 

Let us now state our main results, which are split into four theorems, the two first ones being 
two sub-cases of the Gumbel case, the third one being the Weibull case and the last one the Frechet 
case. In order to state these results, in particular to build some of the normalizing functions g, we 
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need to introduce the following function: for II a measure on integrating x i— > 1 A x and q > 0, 
we denote by yn,<7 the inverse function of the mapping 

tM * J oo (l-e-^)n(dx) + g' (10) 
this inverse function being well-defined on [0, oo) if q > and [(Jq°° xll(dx)) _1 , oo) otherwise. 

Theorem 1.2. Assume that — £ is a subordinator with killing rate q > 0, no drift and a Levy 
measure H such that 

liminf J-^ >0. (11) 

Then I € MDAQumbeb tp = oo and 

m<mi^. | t<T \ — 



t J t— >oo 

Reciprocally, if I £ MDAGumbcl and tp = oo, then necessarily — £ is a subordinator without drift. 

Note that when IT is infinite, the normalizing function t h- >■ t/(pu,q(t) converges to as t — > oo 
and is asymptotically independent of q. But when II is finite, it converges towards 1/(11(0, oo) + </). 
In those cases, the above theorem states that Pi (Xt € • | t < To) has a non-trivial limit. More 
precisely, we have the following equivalence. 

Corollary 1.3. The conditional laws Pi (Xt € • | t < To) have a non-trivial limit as t — )• oo 
z/ and onZy if — £ is a subordinator with no drift and a finite Levy measure. This limit is then 
[ij o (id/ (11(0, oo) + a))" 1 , where II denotes the Levy measure of and q its killing rate. 

The cases were — £ is a subordinator without drift were partly studied in [22] under the as- 
sumptions that q = and the tail distribution of II is regularly varying at with some index in 
[0, 1). The above assumption (jlip on II is clearly far more general than regular variation. When 
(jlip holds the function x t— > J Q X n(y)dy is then said to be of positive increase near since this 
assumption is equivalent to the condition 

J n 2x Ti(y)dy 
liminf " — > 1. 

*-*> j Q x U(y)dy 

We refer to chapter III in Bertoin's book [3] for further equivalent ways of expressing this positive 
increase condition (in particular, see Exercise III. 7 therein). 

In [22], Yaglom limits of non-increasing pssMp were used to describe the asymptotic behavior of 
solutions to a class of fragmentation equations. Theorem 11.21 thus enables us to enlarge significantly 
this class of equations. 

We now turn to the cases of subordinators with a strictly positive drift. The family to which 
belong the Yaglom limit differs significantly according whether the Levy measure is finite or not 
(however, in some sense, the infinite case is the limit of the finite cases as 11(0, oo) — > oo, since 
the exponential distribution is the limit of normalized Beta distribution as considered below). We 
start with the infinite Levy measure cases. Note that the normalizing function in front of X does 
not depend on the killing rate q. 

Theorem 1.4. Assume that — £ is a subordinator with killing rate q > 0, a drift d > and an 
infinite Levy measure II satisfying 

< liminf ; / — < km sup „__; / < 1. (12) 
*->o f Q x U(u)du~ J II(u)du 
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Then tp = 1/d, I 6 MDAcnmbei end 

1 — dt J 1 u / 



Reciprocally, iftp < oo and I G MDAQumbeb then — £ is a subordinator with a strictly positive drift 
and an infinite Levy measure. 

The assumption (|12p on the Levy measure LT is slightly more restrictive than the assumption 
(jlip in the 0-drift case. However it still contains a large class of Levy measures, including those 
with a tail that varies regularly at with some index in (0, 1) (note that we had to exclude the 
index 0). 

In the case of a subordinator with a strictly positive drift and a finite Levy measure, the limiting 
distribution of the re-scaled self-similar process conditioned on non-extinction is no more related 
to ^f 1 . Let B 7 denote a Beta distribution with density 7(1 — i/) 7 ~ 1 l(o,i)(y), 7 > 0, and consider 
the factorization in terms of / of this Beta distribution, that is, for i? 7 independent of /, 

i? 7 / L = w B 7 . (13) 

In Section H] we provide necessary and sufficient conditions in terms of the subordinator — £ for the 
existence of this factor i? 7 . Its distribution is then supported by [0, d] and is uniquely determined 
by its entire moments 

EK]= jjm, „>i, 

i=l ' 

where as before (ft denotes the Laplace exponent of the subordinator — £. We denote this factor 
distribution by 

Theorem 1.5. The process — £ is a subordinator with a strictly positive drift and a finite Levy 
measure if and only if I 6 MDA\y e ibull • I n such cases, 

I t <T ) ^( B (nCo,oo) +9 )/ d ) ; 

d - /Hi 

where as usual q denotes the killing rate, d the drift and LT the Levy measure. 

Remark. This Yaglom limit does not belong to the set of QS distributions of the process X. 

We stress that the proofs of Theorems 11.21 11.41 and 11.51 are based on new results obtained in 
Sections [3] and H] below, some of them being extensions of known results, on the behavior of the tail 
distribution P(J > t) as t — > tp, which are established by using recursively the integral equation on 
the density of / settled by Carmona et al. [13] and extended by [30] . This equation is recalled in 
formula (|16|) , Section 12.11 We believe these results are interesting in themselves and could stimulate 
the study of exponential functionals of subordinators in extreme value theory. 

Last, we focus on the non- monotone cases. We observe that these include the case where the 
underlying Levy process is a killed subordinator. For that end we will need factorizations of the 
Pareto distributions. Let P 7 denote a Pareto distribution with density on R + x 1— > 7(1 + x)~ 7 H 
7 > 0, and consider the factorization 

J 7 / L = W P 7 , (14) 

for some random variable J 7 independent of /. As for the Beta distribution, we will provide 
in Section [5] a necessary and sufficient condition for the existence of such a factorization and 
characterize the distribution of the factor J 7 , which is uniquely determined and denoted by /x]- . 
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Theorem 1.6. The process X is not monotone and admits a Yaglom limit if and only if I € 
MDAp^chot , that is if and only if t ^ P(I > t) is regularly varying at oo, say with index —7 with 
7 > 0. In such a case, 

~r ^ ' I t<T ) —^^ ] . (15) 

Moreover, sufficient conditions for t 1— > P(I > t) to vary regularly with index 7 > as t — > 00 are: 

(i) either that £1 is not lattice and satisfies the so-called Cramer's condition: 

E[exp(7fi)] = 1 and E exp(7£i)] < 00; 

(ii) or Q = co a.s., and 



<->oo P(£ > t) ' ' i^oo P(^ > t) 1 J 

and E[e 7 ^ 1 ] < 1. 7/0 < 7 < 1, is furthermore assumed that E[£i] € (— 00, 0). 
Reciprocally, a necessary condition for I G MDA Fr6chct is tfiat E^ 1 ] < 1 and E[e^ +S ^ 1 ] > 1 /or 
all 5 > 0, /or some 7 > 0. 

The two sufficient conditions (i) and (ii) are proved in [38] and [3D], respectively. We point out 
that in the latter paper the condition in (ii) is stated in an equivalent way for the Levy measure of 
the underlying Levy process. This allows us to ensure that the sufficient conditions for the regular 
variation of t h-> P(7 > t) do not seem far from being necessary as they are the conjunction of the 
necessary condition on the moments of £1 and a rather mild condition on the tail behavior of the 
right tail distribution of the underlying Levy process, equivalently of the Levy measure II. 

Note in passing that the necessary condition on the moments of £1 gives some examples of 
pssMp for which there is no Yaglom limit: e.g. if £ = 07 — o~2 where o"i,o~2 are subordinators, the 
Levy measure IIi of 07 being such that exp(7x)IIi(dx) = 00 for all 7 > 0. 

Conversely, if £ is a Levy process with no positive jumps such that either £ < 00 a.s. or ( = 00 
a.s. and lim^oo £ 4 = —00, then £ has exponential moments of all positive orders, see e.g. [3] 
chapter VII, and hence the conditions (i) in Theorem 11.61 are satisfied. Some specific examples of 
this family of pssMp are studied in Section [71 where other examples are developed. 

Contrary to the monotone cases, the normalizing function g in the Frechet cases is independent 
of the non-monotone process X and is increasing. Note that this linear normalization implies that 
the limiting distribution, when it exists, of 

P*(t _1 ^€- I t<T Q ) 

does not depend on x, since t~ 1 Xt under is distributed as t~ 1 xX tx -i under Pi. Note also that 
this limiting distribution cannot be QS since there is no QS distribution in the non-monotone cases. 
However it is a QS distribution for another Markov process, the Ornstein-Uhlenbeck type process 

U t = e"*X et „ 1 , t > 0, 

which is of interest in itself. Indeed, the above result gives necessary and sufficient conditions for 
the convergence 

¥ x (U t e- I t<TV) 

where denotes the first hitting time of for U. Furthermore, it is well-known that such "stan- 
dard" Yaglom convergence of a Markov process implies that the limiting distribution is QS for the 
corresponding Markov process. In fact we will prove that the existence of a factor J 7 in (|14[) is 
equivalent to the existence of a quasi-stationary distribution for U and we will use this Ornstein- 
Uhlenbeck point of view to describe the distribution y}j 7 . The details of these connection and 
characterization are given in Section [6j 
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2 Background and preliminary results 



This section is divided in two main parts. First some background on exponential functionals of 
Levy processes and a proof of Theorem II .11 are given (Subsections I2.ll and 12. 2| ). Then background 
on extreme value theory are recalled (Subsection I2.3P and applied to set up a "key result" relating 
the existence and characterization of a Yaglom limit of a pssMp X to some asymptotic properties 
of the tail of the associated exponential functional I (Subsection [22 



2.1 Background on exponential functionals of Levy processes 

As we have seen in the Introduction most of our study about existence of Yaglom limits for pssMp 
can be resumed to the study of the residual lifetime of exponential functionals of Levy processes, 
that is why in this section we gather some known distributional results about these r.v. We recall 
that £ denotes a real valued Levy process, possibly killed (in — oo) at time C, and 

I=[ exp(£ s )ds. 
J o 

Assume for the moment that — £ is a subordinator with Laplace exponent (ft, and that its killing 
term is q > 0, its drift is d, and its Levy measure is II. We denote by H(x) the right tail of II, 
Il(x) = II(x, oo), x > 0. Carmona, Petit and Yor [13] established that if q = and the mean of — £ 
is finite then / has a density, say k, and it is the unique density that solves the integral equation 

POD 

k(x) = dxk(x) + / U(ln(y/x))k(y)dy, x € (0, 1/d). 

Jx 

Recently in |30j . it has been proved that both hypotheses above can be removed, in the sense that 
when — £ is a subordinator, the exponential functional / has a density, that we still denote by k, 
and that it solves the integral equation 

roo roo 

k(x) = dxk(x)+ U(ln(y/x))k(y)dy + q / k(y)dy x e (0, 1/d). (16) 

Jx J X 

Conversely, if a density on (0, 1 /d) satisfies this equation then it is a density of /. Carmona et al. 
also obtained the following useful identity about the positive moments of / 



E 



X E 



0(A) 



A-l 



A > 0. (17) 



We now leave aside the assumption that — £ is a subordinator. As we have seen in the Introduction, 
Yaglom limits are qualitatively different depending on whether the support of the law of Tq under 
Pi is bounded or unbounded, so it is important to determine under which conditions the support 
of the law of To under Pi is bounded. The purpose of the following result is to provide necessary 
and sufficient conditions for 

t F = sup{t > : Pi (T < t) < 1} < oo. 

Lemma 2.1. Let £ be a possibly killed real valued Levy process and such that L = Jq e^ s ds < oo 
a.s. We have that the law of L has a bounded support if and only if — £ is a subordinator with a 
strictly positive drift, say d. In this case the support of I is given by [0,d -1 ]. 

Proof. Assume a = — £ is a possibly killed subordinator with a drift d > 0. The bounded variation 
of a implies that it can be represented as a t = dt + X]o< s <( ^-s, t < £, where A s = o~ s — <t s _ > 0, 
for s < C a - s -i see e -g- E] chapter III. Using this representation it is easily verified that 

< 1 - 1 

e~ Ua ds < < - a.s. (18) 

d d 
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Now assume that I has a bounded support that is contained in [0, c _1 ], for some c > 0. This 
implies that / has moments of all orders, and hence using the independence and homogeneity of 
the increments of £ it follows that for any t > fixed and all A > 0, 



oo > E 



> E 



o 

AO 



e^ s ds 
C-t 



o 



E[e A?t ,i < C]E 



Therefore, 1 > E[e A «',t < £], VA > 0, and hence P(& > 0,t < Q = 0. In other case we would have 
that 



1 > lim E 

Afoo 



e x H Ut>0} ,t<C 



E 



lime^l^^o},* < C 

Atoo 



OO, 



which is a contradiction. It follows that P(£t < 0,t < Q = 1, for all £ > 0, which is equivalent to 
say that — £ is a subordinator. 

Let us assume that the drift is equal to zero and establish a contradiction. By Carmona et al.'s 
moment formula (|17p and an elementary upper bound it follows that 



n 

i=i 



i 



E[/ n ] < — , n > 1. 



Taking logarithms in both sides of the inequality and using that for A > 0, = j+ J °° e Xx U(x)d 
with g the killing rate, we get that 

1 n 
ln(c)<-5> 



e~ ix (n(x) + g) dx 



ra > 0. 



Hence using that j + J °° e Xx Il(x)dx > and that the right most term in the latter display is 

A— >oo 

a Cesaro mean, we get by making n — > oo, that 

m(c) = — oo, 

which is a contradiction to the assumption that c > 0. It is therefore impossible to have a zero 
drift whenever the support of / is bounded. Let us now check that the support is given by [0, c£ -1 ]. 
The equation ([18]) implies that the support of the law of /, say [0,c _1 ], is contained in [0,<i -1 ]. 



0(A) 



Repeating the latter argument with — 
therefore d~ l < c" 1 , from where the result follows 



x+ d + Io 



50 e Xx H(x)dx, we infer that lnc < hid, and 

□ 



To finish this subsection we prove the key identity Q. 

Lemma 2.2. Let X be a pssMp that hits continuously, £ be the underlying real valued Levy 
process, and I := Jfj e^ a ds. For < t < tp, 

P(J - t G dy\t <I)= Pi(X t /€ dy\t < T ) 
where I has the same law as I and is independent of (X s , s < t). 
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Proof. Let r(t) = inf{s > : f Q s e^ u du > t}, and observe that r(t) is a stopping time. Conditionally 
on the event {t < 1} = {r(t) < (}, we have 

rC rr{t) ( K-T{t) \ 

I = J e^ds = J e ?s ds + e^W ( J e (^ { t)+u~^(t)) du j = t + e^U, 

where I is independent of (£ s ,s < r(t)) and has the same law as / because the process ((, T ^ +S — 
Cr(i)iO < •? < C — T (0) has the same law as £ and is independent of (£, s ,s < r(t)), by the strong 
Markov property of £. We conclude with Lamperti's identity that 

pr t ,Pi) L = w P). 

□ 

2.2 Quasi-stationary distributions 

The main purpose of this subsection is to prove Theorem ll.il this will be a direct consequence of 
Lemmas 12.31 and !2.4l below. The existence of quasi-stationary distribution for pssMp with decreasing 
paths has been incidentally studied by Bertoin and Yor. An extension of their results is summarized 
in the following Lemma. 

Lemma 2.3 (The case of non- increasing paths. Bertoin and Yor [5]). Assume that X has non- 
increasing paths and hits zero in a finite time, that is the Levy processes associated to X via 
Lamperti's transformation, say £, is such that — £ is a subordinator. Then for 9 > 0, 

(i) there exists a QS-law for X, 

[ fi e (dx) F x (X t e dy\t < T ) = fi e (dy), y>0, 

JR* + 

and such that 

[ ii e (dx)F x (t<T ) = e~ e \ i>0; 

(ii) fie is characterized by its entire moments ((j) still denotes the Laplace exponent of £): 

„ n 

/ x n Mdx) = 0- n Tl>(i), n>l. 



i=l 



(hi) if Re ~ /J-e and Rq is independent of £, then 

-C 



Re [ e^ds L =e/e. 
Jo 



Proof. We assume hrst that 9 = 1. Bertoin and Yor [5] and Berg [2] established the existence of a 
probability measure [i\ such that the claims in (ii) and (iii) hold. The identity for t > 0, 

/ /ii (dx) F x (X t edy,t<T ) = e~ Vi (dy) , y > 0, 

Jr^_ 

has been established by Bertoin and Yor [5] in their Proposition 4 (including the cases when the 
subordinator — £ is killed - this is not clearly specified in their statement, but is clear from their 
"second proof of Proposition 4" ) . Integrating the constant function 1 we get that 

fi 1 (dx)F x {t<T ) = e- t , t>0. 
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From where the identity (i) when 9 = 1. We should now prove the results for any 9 > 0. Let 
fie(dy) := jUi o (# _1 id) _1 (dy). It is plain that the claims in (ii) and (hi) hold. Let / : — > be 
any measurable and bounded test function. To prove the claim (i) observe that by the self-similarity 
of X and the fact that the results are true for \i\ we have the identities 

/ fi e (dy)E v [f(Xt)\t<T Q ]= f ^{dy)E y/e [f{X t )\t<T Q \ 

JR* + JR* + 

= f Ml (dj/) E y [f {9- 1 X te ) \t6 < T ] (19) 

Jr* 

t x 1 (dy)f(9- 1 y) = [ no(dy)f(y). 

JR* + 

□ 

Knowing the results by Bertoin and Yor it is natural to ask if there are other pssMp, which 
have non-monotone paths, that admit a quasi-stationary distribution and if we can obtain a similar 
factorization of the exponential law. The answer to this question is No. 

Lemma 2.4. X admits a quasi-stationary law if and only if X has non-increasing paths. 

Proof. If \i is a QS-law for X, the simple Markov property implies that there exists an index 9 > 
such that 

f n(dx)F x {t<T ) = e~ et , t>0. 
Jr* + 

The self-similarity of X implies 

e~ et = f fM(dx)F x (t<T )= f f i(dx)F 1 (t<xT ). 

JR* + JR+ 

Recall then that (T 0) Pi) L = w (I, P), where I = J C e^ds. Then if R ~ fj, and R is independent of /, 
we have that 

RI L = e/9. 

This identity and the independence imply that / has moments of all orders, thus arguing as in the 
proof of Lemma 12. II we prove that necessarily — £ is a subordinator. This proves the necessity. The 
sufficiency follows from Lemma 12.31 □ 

We have the following corollary to the proof of the latter result. 

Corollary 2.5. Let £ be a Levy process such that I := e^ s ds < oo a.s. There exists an indepen- 
dent r.v. R such that 

-r~. T Law 

RI = e, 

if and only if is a possibly killed subordinator. 



2.3 Some background on extreme value theory 

For details we refer to the books of de Haan and Ferreira [TE] and Resnick Extremal laws are 
the possible limits in distribution of sequences of the form (max{Zi, Z n } — a n ) /b n where the 
random variables Z^i > 1 are i.i.d. and the sequences (a n ,n > 1) and (b n ,n > 1) deterministic, 
with b n > 0, n > 1. When such a limit holds, we say that the (distribution of the) random variable 
Z\ is in the domain of attraction of an extremal law. There are three possible types of limit 
distributions (we say that a r.v. Z is of the same type as Z' if Z aZ' + b for some deterministics 
constants a, b): 
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- Gumbel distribution : ¥(Z < x) = exp(— exp(— x)), x £ ]R 



- Weibull distribution : F(Z < x) = exp(-(-x) 7 ), x < (7 > 0) 



- Frechet distribution : P(Z < x) = exp(— x 7 ), x > (7 > 0). 



For any random variable Z € M, we denote by 



t z = inf{t > : P(Z > t) = 0}. 



The following theorem gathers standard results in extreme values theory. See in particular the 
papers by Balkema and de Haan PQ, Pickands [35] and Theorems 1.1.2, 1.1.3, 1.1.6 and 1.2.5 of the 
book of de Haan and Ferreira |16j . 

Theorem 2.6. For Z a real random variable, there exists a positive function g and a non- degenerate 
probability measure v such that the limit 

P (Z > xgit) + t) 

lim — — = u(x, 00) holds for all x which is a continuity point of z/(-,oo), (20) 

t—*tz ™(Z > t) 

if and only if Z is in the domain of attraction of an extremal law. Moreover, 

(i) if Z is in the domain of attraction of a Gumbel law, then A20\) holds with g{t) = F(Z > 
u)du/F(Z > t), t < tz and u(x, 00) = exp(— x), x > 0. 

(ii) Z is in the domain of attraction of a Weibull law with parameter 7 > if and only iftz<oo 
and x 1 — y ¥(Z > tz — 1/^) is regularly varying at 00 with index —7. Then A20\) holds with 
g(t) = tz — t and v(x, 00) = (1 — x) 7 , x € (0, 1). 

(iii) Z is in the domain of attraction of a Frechet law with parameter 7 > if and only if tz = 00 
and x I—)- ¥(Z > x) is regularly varying at 00 with index —7. Then \20\) holds with g(t) = t 
and u(x, 00) = (1 + x)~ 7 , x > 0. 

Contrary to the Weibull and Frechet cases, there is no simple condition that characterizes 
random variables in the domain of attraction of a Gumbel distribution. For our purpose, we will 
use the so-called von Mises' condition recalled in the forthcoming Lemma 13.31 Other versions of 
von Mises'conditions are also available (see e.g. [36] )• We also mention that the papers of Balkema 
and de Haan [T] and Geluk |2Uj give conditions on conditional moments of the random variables 
(Z — t) + for Z to be in the domain of attraction of a Gumbel distribution, but we will not use this 
approach later. 

2.4 A key result 

Recall the notation X,£,I from the Introduction, as well as notation e, B 7 and P~, with 7 > 0, 
for an exponential r.v. with parameter 1, a Beta r.v. with density 7(1 — x) 7_1 ,x € (0, 1), and a 
Pareto r.v. with density 7(1 + x) _7_1 ,x > 0. The following result links the Yaglom limits of X 
we are interested in to the asymptotic behavior of the residual life-time of /. It is available for all 
Levy processes £ such that I < 00 a.s. 

Proposition 2.7. There exists a function g : M + — > and a non-trivial probability measure fi 
such that 



if and only if I is in the domain of attraction of an extreme distribution. We then distinguish three 
cases : 




(21) 
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(i) I G MDA Gumbel if and only if 



• (I > xg{t) + t) /P(J > t) > exp(-x), x > 0, 

t— >oo 



for some function g. Then we have that 



Xt €■ I t <r )— 



5 (t) 7 

where Hj is the unique distribution such that for R ~ fij independent of I, RJ L = e. 
(ii) J G MDAwcibuii a^d onZy i/iene exists 7 > 0, suc/i i/iat 

P (/ > x(t F ~t)+t) /P(I > t) ► (1 - x) 7 , x G (0, 1). 

£ — >t p 



Then we have that 



-G- t < i J > Hj 



where [i^ 1 ^ is the unique distribution such that for i? 7 ~ independent of I, R^I L = B 7 . 

(iii) I G MDAprechet if and only if there exists a 7 > such that 

P (I > (x + l)t) /P(J > t) > (1 + x) -7 , .x > 0. 

t— >oo 

Then we have that 

'^e- I t<r )-— „?H 

r y t-¥oo 

where fi^f 1 ^ is the unique distribution such that for J 7 ~ l^f 1 ^ independent of I, J^I ^==" P~. 

This result is an obvious consequence of Theorem 12. 6\ combined with the factorization obtained 
in Lemma 12.21 and the forthcoming Lemmas 12.81 and 12.91 Note that it implies the existence of the 
factor distributions ^\ or ^f 1 ^ when, respectively, I is in the domain of attraction of a 

Gumbel, Weibull or Frechet distribution. 

Lemma 2.8. If I is in the domain of attraction of an extreme distribution, there exists some e > 
such that K[I a ] < 00 for all a G (— e,e). 

Proof First note that a random variable in the domain of attraction of an extreme distribution 
necessarily has some strictly positive moments, i.e. E[/ e ] < 00 for some e > 0. This is obvious 
when / G MDAweitmii since its support is bounded. It is also obvious for / G MDApv^cheti since 
t 1— > P(7 > t) is then regularly varying at 00 with a strictly negative index. Last, it is well-known 
that a random variable in the domain of attraction of a Gumbel distribution has positive moments 
of all orders, see e.g. Proposition 1.10 of [36] , 

Under the assumption that the lifetime of the underlying Levy process is infinite, in Lemma 3 
in [30] it has been proved that if E[I 7 ] < 00, for some 7 > 0, then necessarily E[/ 7_1 ] < 00. But the 
result is true in general and the proof in [30] can be easily extended. In particular Ef/ 6 ^ 1 ] < 00, 
which leads to the claimed result, taking e G (0, 1/2) if necessary. □ 
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Lemma 2.9. Let U, Z,Yt,t > be strictly positive random variables, with Z independent of (y , t > 
0) and such that K[Z a ] < oo for all a G (—e,e) for some e > 0. Then we have the implication, 



Y t Z U=>3a r.v. Y such that Y t y 

t— »oo t— >oo 



where the distribution of Y is uniquely determined by the equation 

YZ L = U 

for Y independent of Z. 

Proof. If K[Z a ] < oo for a G (—e,s) for some e > 0, then z G C h-> Efe"' 11 ^] is analytic on 
{z £ C : Im(z) G (— e, e)}. In particular, .4. := {i € R : E[e* tln ^] = 0} has no accumulation point 

and there exists a neighborhood of which it does not intersect. From the convergence YtZ U, 
we get that 

E[e mlnYt ] > J ■ , J Va G R\A 

Since the limiting function is well-defined and continuous on a neighborhood of 0, we obtain, exactly 
as in the proof of Paul Levy's continuity Theorem, that (lny,/; > 0) is tight. Let lny and hxY' 
be two possible limits in distribution. Then E[e ialny ] = E[e ialny '] for all a G R\A Since A has 
no accumulation point and the characteristic functions are continuous on K, there are identical. 
Hence lny is distributed as lny', and lny converge in distribution. We denote by lny the limiting 
random variable. Necessarily, there exists then a version of Y independent of Z and such that YZ 
is distributed as U. Moreover, if Y' is a strictly positive random variable independent of Z and 
such that Y'Z L ^ U, then the characteristic functions of lny and lny' coincide on R\A Hence 
they are identical and lny' is distributed as lny. □ 

3 Yaglom limits: Gumbel cases 

The goal of the present section is to prove Theorem 11.21 Corollary 11.31 and Theorem 11.41 In view of 
Proposition 12.71 our goal is to characterize the self-similar Markov processes whose extinction time 
belongs to the domain of attraction of a Gumbel distribution, and then determine the normalizing 
function g. Necessary conditions are easy to settle: it is well-known that if / G MDAQumbei, then 
it possesses positive moments of all orders, which implies that — £ is a subordinator, as observed in 
Section 12.11 We will moreover see in the next section that when — £ is a subordinator with a drift 
d > and a finite Levy measure, I is in the domain of attraction of a Weibull distribution. From 
this we conclude that necessary conditions for I G MDAcumbei are: 

• either that — £ is a subordinator without drift (then = oo) 

• or that — £ is a subordinator with a strictly positive drift and an infinite Levy measure (then 
tF < oo). 

Reciprocally, the two following propositions give sufficient conditions for I G MDAQumbei- We recall 
that fn,q is defined from a Levy measure n and a real number q > via (jlOp , 

Proposition 3.1. Assume that — £ is a subordinator with killing rate q > 0, no drift and a Levy 
measure H satisfying (TI\). Then, 

I > t + 



> expf— x), \fx > 0. 

P(J > t) t^oo FV ; ' 
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Proposition 3.2. Assume that — £ is a subordinator with killing rate q > 0, drift d > and an 
infinite Levy measure II satisfying I7g|) , Then, 



I >t + iyiiu ( -r-j 

> exp(— x), Vx > 0. 



> t) 1 



More precisely, together with Proposition 12.71 these two propositions and the above discussion 
prove Theorems 11.21 and 11.41 It remains therefore to prove these propositions, and for this we will 
use the so-called von Mises' condition, which is reminded in the following lemma. 

Lemma 3.3 (Resnick [36], Prop. 1.17). Let U be a non-negative random variable with density f 
such that the von Mises' condition is satisfied 

f(t) f°°¥(U>u)du t , 

\ku>t)Y -^ s ^>0:P(C/>r)>0}. 

Then, as t -> sup{r > : P(U > r) > 0}, 

P (U > t + x / t °° P(f/ > u)du/P(^ > *)) 



> exp(— x), x > 0, 



P(C7 > t) t^oo 

and a similar convergence holds when replacing J^°P(U > u)du/P(U > t) by any asymptotically 
equivalent function. 

We recall that under the assumptions of Propositions 13.11 and 13. 2| the random variable / pos- 
sesses a density k which satisfies the equation f)16|) . By using recursively a variant of this equa- 
tion, we will obtain the same estimates in terms of <pn,q f° r both quotients k(t)/¥(I > t) and 
P(/ > t)/ A oc P(7 > u)du as t — > tp. This will prove von Mises' condition for 7, hence Propositions 
l3.1l and l3T2l Besides, by integrating k(t)/¥(L > t), we will also get an asymptotic estimate in terms 
of (fu,q of — lnP(I > t), which extends previous results on this topic. 

For the proofs of these estimates, we separate the drift-free cases (Section 13. 1|) from the strictly 
positive drift cases (Section l3.2|) . In all cases, <p denotes the Laplace exponent of the subordinator 

-e 

Last, Corollary 11.31 is proved in Section | 



3.1 Drift-free case: proof of Proposition I3TT1 

In this section it is assumed that — £ is a subordinator with killing rate q > 0, no drift and a 
Levy measure denoted by II. The random variable / then possesses a density k which satisfies the 
equation (|16p by Carmona, Petit and Yor with d = 0. Integrating by parts, this becomes 



k{x) = I U(ln(u/x))k(u)du + qP(I > x) 

J X 

roo 

= I > x)-¥(I > xe v ))U(dv)+q¥(I > x), (22) 

Jo 

which, setting /(x) := — ln(P(J > x)), leads to the equation 

/'(*) = = J™ (i - ex P(/(^) - n(d^) + q. (23) 
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Proposition 3.4. Assume moreover that H satisfies All]) - Then, 



\ X — »OG 

> X) X 



/'(*) : 
and consequently 

f(x) = - ln(P(7 > x)) f ^^dn, 

Ja U 

a being any real number such that (fu,q is defined on [a, oo). 

In particular, when is regularly varying at with index —f3 € ( — 1, 0] the assumption pip 
on II is satisfied and standard Abelian-Tauberian theorems [9] imply that <pn,q varies regularly at 
oo with index 1/(1 — /?) and then f(x) ~ (1 — /3)(^n,(j(^)- We therefore recover the second author 
estimates for — lnP(I > x) [37], with a different approach: 

Corollary 3.5. When Ii{x) is regularly varying at with index —(3 £ (—1,0], then 

- ln(P(7 > x)) -^oo (1 - P)<pn, q {x). 

Remark. The assumption (fTT|) on II is only needed to get that limsup^^ (f'(x)x/ipn,q(x)) < 1. 
Moreover, we will see in the proof below that when II is infinite, without any further assumptions, 

f'(x) > ipYi,q(x)/x for all x large enough. 
Next, we set up a similar estimate for ¥(I > x)j f^°F(I > u)du. In that aim, set 

POO 

~Fi(x) = I P(J> u)du, and g(x) = -lnF/(a;). 

J X 

Note that Fj(x) < oo since E[7] < oo. Then integrate equation (|22p to get, together with Fubini's 
theorem, 

P(J > x) = J (jl(x) - e- v T I {xe v )^j H(du) + qFi(x) 

and therefore 

w(t > r ) r°° 

9'^ = r°°wr ^ \a = & ~ exp (~ V + 9 ^ ~ g ( xeV ^ U ( dv ) + 9, x > 0. (24) 
J x P(J > u)du J 

Although this equation is not exactly the same as that satisfied by /, we will obtain in a similar 
manner the following estimate. 

Proposition 3.6. Assume tll\) . Then, 

, _ P(/ > X) m,q{x) 

J P(J > u)du x 
Proposition 13.11 is then proved by combining Propositions 13.41 and 13.61 together with Lemma 

E3 
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3.1.1 Proof of Proposition 13.41 

A first easy observation, is that x *— > ¥(I > x) is rapidly varying at oo in the drift-free case (without 
any assumption on II, except that 11(0, oo) > 0). This means that for any c > 1, 

P(J > cx) 

—— r > as x — >• oo. 

P(I > x) 

This is easy to see. Note first that for c > 1 and all x > 0, 

P(/>cx) < f™F(I>u)du /~P(/ > u)Au 1 



> x) ~ x(c - 1)P(J > x) ~ x(c - 1)P(J > x) x(c-l)g'(x) 
But d2U) implies that g'(x) > 0(1) > for all x. Hence, 

P(7 > cx) < 1 



P(J > x) ~ x(c- 1)^(1) 
and letting x — > oo, rapid variation is proved. 

(1) II is finite. The proof of Proposition 13.41 is easy when II is finite. On the one hand, we 
immediately get from the definition of tpn,q, that tpn,q(x) ~oo x(II(0,oo) + q). On the other hand, 
the rapid variation of x i— > P(7 > x) at oo implies that exp(/(x) — f(xe v )) — > as x — > oo for all 
v > 0. Plugging this in ([23]) and using Fatou's lemma, we get that 11(0, oo) + q < liminf r _ 5 . 00 f'(x). 
That f'{x) < II(0,oo) + q for all x > is obvious. Hence f'(x) ~ n(0, oo) + q ~ cpu,q(x)/x as 
x — > oo. 

(2) n is infinite. We will need the following remark. Let (0'(O+)) _1 := lim^o x/<p(x) < oo. Since 
x i-> x/4>(x) is continuous strictly increasing from (0, oo) to (((^'(O-l-)) -1 , oo), the function yn,g is 
well-defined on ((</!/(0-|-)) -1 , oo) and for x > ((//(0+))"" 1 , ipu,q(x)/x is then the unique non-zero 
solution to the equation (in u) 

u = 4>{xu). (25) 

Note that this is true for all (finite, infinite) Levy measures n. But the following lemma is only 
true for infinite Levy measures. As above, it relies on the rapid variation at oo of the tail P(I > x) 
and on Fatou's lemma, and so we omit the details. 

Lemma 3.7. As x — > oo, f'(x) — > oo. 

Now, we construct by induction a sequence of non-negative continuous functions h n defined on 
[0, oo) by ho(x) = l,Vx > and 



oo 



exp ^— J h n (u)duj j n(dt>) + q, x > 0. 



h n+ i(x) = J ^ 1 - exp ( - / h n (u)du ) ) U(dv) + q, x > 0. (26) 
Then, 

• for each n > 0, h n is non-decreasing on (0, oo): this is proved by induction on n, using that 
when h n is non-decreasing, so is x i-4- J^ e h n (u)du; 

• for each n > 0, and all x > 0, h n+ i{x) > (f>(xh n (x)). To see this we use that h n is non- 
decreasing, hence x h->- J^ e h n (u)du > h n (x)x(e v — 1) > xh n (x)v. The conclusion follows 
from the definition of d>\ 
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for all x > xq for some xq large enough, the sequence (h n (x),n > 0) is non-decreasing and 
its limit, denoted by h OQ (x), satisfies h OQ (x) > ipu,q(x)/x. Indeed: since hi(x) > <f>{x) and 
since 4>(x) has an infinite limit as x — >• oo, we have that h\{x) > Hq{x) for all x > xq for 
some xq < oo, that is moreover assumed larger than ((/>' (0+))~ 1 . By an obvious induction, we 
then have that /i n +i(x) > h n {x) for all x > xq and all n > 0. Hence the existence of h 00 (x). 
Last, if hoo(x) = oo, it is obviously larger than (pu,q(x) / x . Otherwise, use the fact that <f> is 
continuous and h n+ \{x) > (b(xh n {x)) for all n, to get that hoo(x) > 4>(xh 00 (x)) for x > xq. 
Note then that h oc (x) > ho(x) = 1 (hence is non null), to get that xh 00 (x)(4>(xh 00 (x)))~ 1 > x, 
hence that xh QO (x) > (pn,q(x) on [xo,oo) } using for this conclusion the few lines preceding 
equation (f25l) . 



Now, by Lemma 13.7} there exists x\ such that f'{x) > 1 = ho(x) on [x\,oo). By induction on 
n, since /' satisfies equation ([23|) . we easily get that f'(x) > h n (x) on [xi,oo) for all n > 1. Hence 



We conclude that 



We will next prove 



f'{x) > h OD (x) > ifu, q (x)/x on [max(x ,xi),oo) 
liminf (f'(x)x/tp u , q (x)) > 1. 



limsup (f'(x)x/ipn, q (x)) < 1, 

a;— >oo 

which turns out to be more tricky. It is a consequence of the two following lemmas. 
Lemma 3.8. f'{x) = hoo(x) for x € [xq,oo). 

Proof. It relies on the following observation: / is the unique differentiable function defined on 
[0, oo) which is non-decreasing, null at 0, converging to oo as x — > oo and satisfying (|23p . Indeed, 
let h be a function satisfying to all these properties. Then x € [0, oo) — > 1 — exp(— h(x)) is the 
cumulative distribution function of some random variable possessing a density satisfying equation 
(|22p . But, as already mentioned, this equation has a unique solution which is a density. Hence 
h = f. 

Now, consider an extension of defined on [0, oo): for this set ho(x) = h^x) on [xo, oo) and 
ho(x) = on [0, xo). Then construct a sequence (h n ,n > 0) by induction, constructing h n+ \ from 
h n as h n+ \ was constructed from h n via (f26|) . Clearly, h n (x) = h OQ (x) on [xq, oo) for all n. Besides, 
for all x € [0, xo), hi(x) > = ho(x) and of course, this also holds for x > xq. Hence h\(x) > ho(x) 
on [0, oo) and this implies, by an obvious induction, that (h n (x),n > 0) is an increasing sequence 
for all x > 0. Call h OQ (x) its limit for all x > and let be its primitive on [0, oo) nul at 0. 
Clearly: is non-decreasing, differentiable, converges to oo as x — > oo and satisfies (f23|) . Hence 
Hqo and /ioo are uniquely determined and hoo = f ■ □ 

Lemma 3.9. Assume 07]), Then limsup x ^. 00 (h OQ (x)x/ipu,q(x)) < 1. 

Proof. A key point of the proof is that (|11|) implies the existence of reals numbers < (3 < 1 and 
x\ > such that for x > x\, 

<p(ax) < a p <f>(x), for all a > 1. (27) 

To see this, the first step is the equivalence between (jlip and 

t^'(t) 
lim sup — rr- < 1 

t^oo q>(t) 

(see e.g. Exercise 7 of - and note that it is still true when q > 0). Hence the existence of some 
(3 € (0, 1) such that for all a > 1, and all x large enough, 
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which indeed gives <p(ax) < a^(f>(x). Take then /3' € (j3, 1). This implies in particular that 

(j){x) < x 13 ' for all x>xi, (28) 

taking x\ larger if necessary. 

We will need the following fact. Consider e € (0, 1) and use that for x large enough and all n, 

4>{x) < K +1 (x) < J ( 1 - exp ( - J h n (u)duj J n(du) + n(ln(l + e)) + q, 

together with the fact that Il(ln(l +e)) < e4>(x) < eh n+ i(x) for large x (since II is infinite) and all 
n, to get that 

/in+i(x) < (l-e)" 1 (jf ( + } (l-exp(-J h n (u)d u y^ U(dv) + 

< (1 - e)" 1 ^ ( +£) - exp (-x ln(i £ +g) ^ ra (x(l + e))) ) n(dt;) + ^ 

where we have used for the second inequality that h n (u) < h n (xe v ) < h n (x(l + s)) for u < xe v and 
v < ln(l + e), since h n is non-decreasing, and also that v~ 1 (e v — 1) is increasing on (0,ln(l + e)). 
Then for x large enough (depending on e), say x > x £ (take x e > max(xi,l) for simplicity), one 
has 

-xh n (x(l +e))) , for all n > 0. (29) 



K+i(x) < (l-e)" 1 ^ 



>(l + £) 

Now, consider the sequence of functions defined inductively by 4>o(x) = x^'/C 1 -/ 3 ') anc l for n > 1, 

4> n (x) = 4)(x(j) n -i(x)), x > 0. 

Using (j28|) . we see that </>i(x) < x^+^'/C 1 -/ 3 '))/ 3 ' = <^o(x), Vx > xi. Thus, for those x, (4> n (x),n > 1) 
is non- increasing (and non-negative). Its limit 4>oa(x) satisfies (froo(x) = <j)(x(j) 00 (x)). According to 
the lines above ([25]) . this implies that either <f)oo{x) = or </>oo(x) = ipn,q(x)/x. 

Besides, by induction, using (f27|) and that the functions 4> n are all non-decreasing on [0, oo), we 
obtain, 

M{l + e)x) < (l + e)^-i^ + ^> n (x), x>0,n>0, (30) 

with 7 = P'/(l - /3'). Next, note that h (x) = 1 < x^7(i-/9') for x > 1. Hence, using ((29)) together 
with (j27|) and (|30p . we get by induction on n that for all x > x e and all n > 

Then, for x > max(xo,x e ), letting n — > oo, 

) (mTTT)) d + O^'^feW. 

Among other things, this implies that (j)oo( x ) > (we remind that hoo(x) > h-o{ x ) = 1) hence 
0oo(x) = ipYia{x)/x, and 

/ioo(x)x / 1 \£i>°^ / e \£i>i/3 l 
limsup y < — r (l + e)^'> l(J ^ , 



^n,g(a;) V 1-6 / Vln(l+. 
and the upper bound converges to 1 as e tends to 0. □ 
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3.1.2 Proof of Proposition |3"U1 

(1) II is finite. We have seen in the previous section that when II is finite, k(x) ~oo (11(0, oo) + 
q)f(I > x) as x — > oo. Integrating this equivalence, we get that k(u)du ~oo (11(0, oo) + 
q) Jf°P(7 > u)du as x — > oo. Hence </(x) ~oo (pn. q {x)/x. 



(2) II is infinite. The main lines of the proof are similar to those of the proof of Proposition 13.41 
when II is infinite. We only detail the main differences. First, for A > 0, we know by Lemma 13.71 
that k(x) > A¥(I > x) for x > xa- Integrating this inequality , we get that g'(x) > A for x > xa, 
hence g'{x) — ► oo as x — > oo. 

Introduce next the sequence of functions h n ,n > defined recursively by 

h n +i(x) = J ^1 - exp f—v - J h n (u)du^j^j U(dv) + q, x>0, (31) 

starting from ho constant, equal to 1. Noticing that h n+ i(x) > + xh n (x)) > (f>(xh n (x)) for 
x > 0, n > 0, we obtain, exactly as in the proof of Proposition 13.41 using ()24p . that 

g'(x) > hoo(x) > ipu,q(x)/x, for large enough x, 

with obvious notations. 

For the converse direction, we need the following lemma. Recall the notation f3' introduced in 



Lemma 3.10. Assume ill]). Then for all x large enough, 

g\x)<2{2xfl^'\ 

Proof. First, note from (|28p that for x large enough, ipu tq (x)/x < x^'^ l ~^'\ Hence, by Proposition 
EH f'{x) < x 13 V(l-^')3/2 and then k(x) < x^'/^^'^(I > x)3/2 for large x. Integrating this 
inequality on [u, 2u] for large u, we get 

o poo 

P(J > u) - P(J > 2u) < -(2n)^'/(i-/3') / p(/ > y )dy. 

2 Ju 

Besides, since u i— > P(7 > n) is of rapid variation as n — ?■ oo, then P(7 > 2u)/P(7 > u) — > as 
u — > oo. Hence the conclusion. □ 

Last, it is easy to prove that 
Lemma 3.11. Assume if77]). TTien lim sup^^ (g'(x)x/(pu,q(x)) < 1. 

Proof. Introduce the sequence of functions h n ,n > defined recursively by ho(x) = 2(2x) /3 ) 
and for n > 1, 

h n+1 (x) = J^ (l-exp(-v-J h n (u)du\) U(dv) + q, x>0. (32) 

Trivially, using Lemma 13.101 and (|24p . we get that for all x large enough such that g'(x) < 
2(2^7(1-/3'), 

^'(x) < h n (x), for all n > 1. 



The proof is then very similar to that of Lemma [3 .91 except that we have to replace inequality (|29p . 
by 

h n+1 (x) < (1 - e)- 2 4> I — -xh n {x(l + e)) ] , for all ?i > 0. (33) 

\ln(l + e) J 
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The reason for this is that (1321) leads first to 



xh n (x(l +£)) + !), for all n > 0, 



>(l + e) 

for all x large enough, independent of n. But from ()27|) . we have that 

0(1 + a:) < (l-e)"V(^), 

for x sufficiently large (depending on e). Since, (e/ln(l + e)) > 1 and h n (x(l+e)) > g' (x(l+e)) > 1 
for all n and all x larger to some xq independent of n, we indeed get (|33[) . 

As in the proof of Lemma f3T9l keeping the notation /3,7 introduced therein, this leads to 

for all x sufficiently large (depending on e but not on n), where (4> n , n > 0) is a non-increasing 
sequence of functions defined recursively from the same induction relation as the one used to 
construct the sequence (4> n , n > 0), except that the initial function is here 4>q = h$. The conclusion 
follows similarly. □ 

3.2 Positive drift cases: proof of Proposition 13.21 

We now assume that the subordinator — £ has a drift d > and an infinite Levy measure II. In 
this case Carmona, Petit and Yor's equation (|16p for k becomes 



(l-dx)k(x) = / U(ln(u/x))k(u)du+qF(I > x) = (P(J > x) - P(J > xe v )) Ii{dv)+qW{I > x), 



o 



for a; € [0,1/d) and k(x) = otherwise. Setting, as in the previous section, f(x) = — ln(P(7 > x)), 
we get that for x £ [0,1/d) 

/•oo 

(l-dx)f'(x)= (l-exp(f(x)-f(xe v )))H(dv) + Q (34) 
J o 

and f{x) = oo on [1/d, oo). Note that for all v > 0, P(I > xe t, )/P(I > x) is equal to for 
x € (e~ v /d, 1/d). Hence by Fatou's lemma, we see that since II is infinite 

n t \ .// \ (1 - dx)k(x) 

1 - da;)/ (x) = > oo, as x -)■ 

F(i > 

In addition to their application to Yaglom's limits of positive self-similar Markov processes, the two 
following results are also interesting by themselves as they give new results on the tail behavior of 
P(J > x) as x tends to the right most extreme of the support of I when this support is bounded. 



Proposition 3.12. Assume that holds. Then, 



f'(x) ~ dipu,o \ TZ~d^) ' cis x \ 1/d, 

and consequently 

f(x)~dj yn,o ( - dtt ' asx^l/d, 

where a G (0, 1/d) is taken large enough so that (pu,o(u/(l — du)) is well-defined on (a, 1/d). 
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Corollary 3.13. //II is regularly varying at with index f3 6 (0, 1), 

Next, consider as in Section [3,11 the function g(x) = — \n(^f^°¥(I > u)du), which is equal to 
co on [l/d,oo), and which, integrating Carmona, Petit and Yor's equation, satisfies on [0, 1/d) the 
equation 

poo 

g'(x)(l - dx) = d+ / (1 -exp(-v + g(x) - g(xe v ))U(dv) + q. (35) 
Jo 

Proposition 3.14. Assume that U2\) holds. Then, 

g'(x) ~ d<p n ,o \ TZ~^) ' as x ^ l / d ' 
Together with Lemma 13.31 this proposition and Proposition 13.121 imply Proposition 13.21 
3.2.1 Proof of Proposition l3?L2l 

Again, the proof is similar to that of Proposition 13.41 First, using that (1 — dx)f'(x) — > oo as 
x 1/d, we easily obtain, mimicking the proof of Proposition 13.41 that for x large enough 



f'(x) > aVn, g (^) ~, t i/ d dm,o (^) 



(recall (pu,q(x) ~oo ^Pu,o(x) when II is infinite). This is done by considering the sequence of functions 
h n defined by ho(x) = 1/(1 — dx) on [0, 1/d), ho(x) = oo on [1/d, oo) and then by induction 



h n+ i(x) = (1 - dx) 1 J ^l-exp^y h n (u)duj^j 



h n (u)duj J U(dv) +q, x£[0,l/d) 

and h n -\-\{x) = oo for x > 1/d. We claim that for x 6 [xo,l/d) for some xo < 1/d large enough, 
the sequence (h n (x),n > 0) is increasing, with finite limit h OQ (x) < f'(x), and moreover, 

K+i(x) > 1 4>(xh n (x)), Vn > 0, 
1 — ax 

where cft(x) = 4>{x) — dx = f^°(l — exp(— xy))LT(dy) + q. This leads to the expected lower bound, 
details are left to the reader. 

The other direction cannot be adapted so easily. However the proof of the following lemma is 
really the same as that of Lemma 13.81 using the uniqueness of a density solution to the Carmona, 
Petit and Yor's equation. 

Lemma 3.15. f'(x) = h^x) for x € [xq, 1/d). 

Hence Proposition 13.121 will be proved once we have checked the following result. 
Lemma 3.16. Assume U2\). Then limsup^^-iy^ (jioo(x)/<fu,g ( i-dx ^ ) — d - 
Proof. For < e < 1, the correct way to split the integral defining h n+ \ from h n is now 

h n+1 (x) < (1 - dx)- 1 i (l-expf-/ h n (u)dujjU(dv) + U(ln((dx)- £ ))+q\. 

(36) 
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We will use that U(ln((dx) £ ) is negligible compared to h n +i(x) as x — >• 1/d. To see this, note from 
the discussion above that for n > 1, 

1 — ( x 

K+i(x) > h 2 (x) > 



1 — dx \ 1 — dx \ 1 — dx 

for x G [xo,l/d). Moreover, the right-hand (strict) inequality in the assumption (|12p implies that 
liminft-^oo tcj) (t)/(j)(t) > (this is an easy exercise), which in turn leads to the existence of a 7 > 
such that 4>{ax) > a 7 0(x) for all a > 1 and all x large enough (this is proved similarly to ([2j 
Hence for all A > 0, for x sufficiently close to 1/d (depending on A) and all n > 1, 

K+i(x) > h 2 (x) > 



1 — dx \ 1 — dx 

Besides, it is a classical result - that can easily be checked - that U(u) < C(p(l/u) for some constant 
C when u is close to 0. Also, 4>{\u) < X(j)(u) for A > 1, by concavity of </>. From all this we deduce 
that for x sufficiently close to 1/d (independently of n, but not on e), (1 — dx)~ 1 H(ln((dx)~ £ ) < 
e(l + e)~ l h n+ i(x), and therefore, from ([36]) . 



l + e f fH(dx)-) 
h n+ i(x) < 1 _ dx 



1 - exp (- j h n {u)Au\ U(dv) + q 



< T^rfi (h n {x(d X y e ))x 



1-dx V y y ' " \ln{(dx)~ £ )) _ 

where we have used in the second inequality that the function v G i— >• (e v — l)/v is increasing. 
Moreover, taking x closer to 1/d if necessary, we have that ((dx)~ e — 1)/ ln((dx)~ £ ) < (l + e), which 
finally leads to 

K+i(x) < £L±fL-0 (h n (x(dx)- £ ))x) (37) 

for some < (3 < 1, and all x G [x £ , 1/d) for some x £ < 1/d, and all n > 1, using that under (fT2j) . 
inequality (j27|) is valid by replacing there the notation by (of course the same remark holds for 
(f28|) ). We also keep from there the notation j3' (see ([28]) ). 

Next set 4>q(x) = d~@ )(1 — dx)™ 1 ^ 1 "^ ) for £ G [0, 1/d) and define recursively, 

0n+l(g) = 1 — d J > (4>ni x ) x ) , X<l/d. (38) 

Using (f28j) for 0, we see that 4>i(x) < <f> (x) for x sufficiently close to 1/d, hence the sequence 
(4> n (x),n > 0) is non-increasing on this neighborhood of 1/d. Its limit (^(x) satisfies 4>oo( x ) = 
(1 — dx)~ l (f) ( y cf) 00 {x)x) , hence is either equal to or x" 1 ipn,q(x / (1 — dx)) (since ipn,q is the inverse 
of x 1-4 x/(j){x)). 
Last, note that 

1 1 

< x 



1 - {dx) 1 - 6 ~ 1 - £ 1 - da: 

on a neighborhood of 1/d, as well as (dx)~ £ < 1 + £, still for x near 1/d. Combining this with ([3 
and (|38p we get, by induction on n, that on a left neighborhood of 1/d (depending on e) 

h M 1 ■ /1 -i- AlS'^/M 



which, letting first n — > 00, then x — > 1/d and then e — >• 0, leads indeed to 

limsup (hoo(x) /<pn,q (x/(l — dx))) < d. 



□ 
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3.2.2 Proof of Proposition ETT31 

Again, the proof holds in a way which is similar to the proofs of Propositions 13. 4[ 13.61 and 13.121 
although some adjustments are necessary. We point out that the drift d in the right-hand side of 
equation ([35]) plays a negligible role. We only detail precisely here the starting (rough) bounds for 
g' which are sufficient to initialize the inductions. 

On the one hand, we easily get from equation (|35p that g'(x)(l — dx) — > oo as x — > 1/d. This 
comes from the fact that for v > 0, g{x) — g(xe v ) = — oo provided x < 1/d is sufficiently close to 
1/d and also from the fact that IT is infinite. This is enough, using ([35]) . to get without difficulty 

that g'(x) > dipu^q f id^fa ) ^ or x l ar S e enough. 

On the other hand, note that Proposition 13.121 and assumption (j!2|) imply that f'{x) < 
(1 — <ix) -1// ( 1- ' 3 ') for some j3' £ (0,1) and all x sufficiently close to 1/d (since (|12p implies an 
inequality of type ([28]) for cj>). Hence, for those xs, 

[•x^xd)^ 1 / 2 poo 

>x)- P(J > x(xdy 1/2 ) = / k(u)du < (1 - {dx) l/2 )~ ll{l - p,) / F(I > u)du. 



Moreover, P(I > x(xd)~ l l 2 ) /F(I > x) — > as x — > 1/d. Indeed, since (1 — dx)k(x)/F{I > x) — » oo 
as x — > 1/d, we get that for all a > 0, aP(J > x) < (1 — dx)k(x) for x sufficiently close to 1/d. 
Then for x large enough, x < 1/d, 

P(J > x(xri)- 1 / 2 ) ff (gd) ~ 1/2 P(J>tt)dtt ^ (xd) " 1/2 (l-^)fc(^)d^ 

P(/ > x) - x((xd)~ 1 /2 _ i)p(/ > x ) ~ ax((xd)-V2 _ y-py > x ) 

rxixd)- 1 / 2 



(l-dx)f* (xa> k(u)du < 1-dx 



(IX 



((xd)-V2 _ i)p(/ >x)~ ax((x<i)- 1 /2 _ i) ' 



Letting x — > 1/d, we get that 



p(/ > xixdy 1 / 2 ) 2d 

llmSU P tc/T- ^ \ ^ — ' 



for all a > 0. Hence P(J > x(xd)" 1 / 2 )/P(J > x) -> 1/d. From all this we deduce that 

</(x) < C(l - (dx))- 1/{1 ~^ 

for x close to 1/d, with C > 1. But this is a sufficiently nice upper bound for g' to use the usual 
schemes to get that 

limsup (g'(x)/ipu, q (x/(l - dx))) < d. 

x-^l/d 

The idea is to construct inductively a sequence of functions {h n ,n > 0) with an induction scheme 
based on equation ([35]), starting from ho(x) = C(l — (dx))~ 1 ^ 1 ~^'\x < 1/d. Then, clearly, 
g'(x) < h n {x) for all n. But also for n > 2 and x sufficiently close to 1/d, 

I _ ( x — 

hJx) > 



1 — dx \ 1 — dx \ 1 — dx 

which can be proved inductively using that ho(x) > (1 — dx)^ 1 . This is sufficient to settle an 
inequality similar to ([37|) for the sequence {h n , n > 1). Last, it is easy to adapt the end of the proof 
of Lemma 13.161 to get the expected upper equivalent function. 
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3.3 Proof of Corollary Q 

Part of this corollary is a consequence of Theorem 11.21 The only remaining thing to prove is 
that the convergence of ¥\(Xt 6 -\t < To) to a non-trivial limit as t — > oo implies that — £ is a 
subordinator with drift and a finite Levy measure. From Proposition 12.71 we already know that 
this convergence implies that I € MDAcumbei with tp = oo, hence — £ is a subordinator with 
drift, and that there exists a constant c > such that for x > 

P( J > cx + t) . . 

— ^tt ; > exp — a;). 

Suppose now that the Levy measure of — £ is infinite and recall from Lemma 13.71 that this implies 
that for all a > 0, k(u) > aP(I > u) for all n large enough. Hence for t large enough (x is fixed), 

P(J > cx + 1) /, ca:+f P(/>n)dn J t cx+t fc(u)dM 1 
P(7 > i) ~ cxF(I > t) ~ acxP(I > t) ~ acx' 

And therefore exp(— x) < (acx) -1 for all a > 0, which is absurd. Hence the Levy measure of — £ is 
finite. 

4 Yaglom limits: Weibull cases and factorizations of Beta distri- 
butions 

In this section, we prove Theorem 11.51 and give a necessary and sufficient condition on the distri- 
bution of the Levy process £ for / to be a factor of the random variable B 7 . 

4.1 Proof of Theorem 11.51 

Recall from Lemma 12. II that tp = 1/d when — £ is a subordinator with drift d > 0. We start with 
the following result. 

Lemma 4.1. When — £ is a subordinator with killing rate q > 0, drift d > and a finite Levy 
measure U, 

1 A k(t) ^n(0,oo) + g 



d ) P(/ > t) t_>i d 
where k denotes the density of I. 

Proof. For t < 1/d, Carmona, Petit and Yor's equation (|16p becomes 

1 t )w-± r(i- r 4P^)wv) + 



d /P(J>t) rfj \ P(J>t) / v ' d 
For all v > 0, P(7 > te 1 ')/P(I > t) is equal to for t sufficiently close to 1/d, hence 

1 t \ k(t) n(0,oo)+g 



d J F(I > t) t-n/d d 
by dominated convergence. □ 
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Proof of Theorem \1.5[ • Assume that — £ is a subordinator with killing rate q > 0, drift d > and 
finite Levy measure II. It is then a standard result of regular variation theory that the convergence 
established in Lemma 14.11 implies that for x € (0, 1), 

F(I > t + (1/d - t)x) . , n(o,oo)+„ 

— s > (1 — a;) d 

P(J>t) 1 1 

and we conclude with Proposition 12.71 (ii) . 

• Reciprocally, assume that I € MDA-wcibuii- Then, / has a bounded support, which implies 
by Lemma 12.11 that — £ is a subordinator with a strictly positive drift d. Moreover, the function 
u — > P(J > 1/d — 1/u) is regularly varying with an index —7, 7 > 0, as u — > 00. When 11(0, 00) is 
infinite, this is incompatible with the fact (1/d — t)k(t)/¥(I > t) — > 00 as t —¥ 1/d, proved at the 
beginning of Section [3.2i Hence 11(0, 00) is finite. □ 

4.2 Factorization of Beta distributions 

In the sequel 7 will denote a Beta random variable with parameters /3,7 > 0, and when (5 = 1 
we will simply denote it by B 7 . 

Proposition 4.2. Let 7 > 0. There exists a random variable i? 7 independent of I such that 

T E> "R 

1 Ity XJ~ 

if and only if is a subordinator with a strictly positive drift d and a finite Levy measure U such 
that LT(0, 00) + q < d*y. In such a case, sup {t > : P(i? 7 > t) > 0} = d and the distribution of Ry 
is characterized by its entire moments which are given by 

E[Ry = f[^L, n>l, (39) 

with (j) the Laplace exponent of — £. 

Proof. • Assume that the factorization R^I L = ¥ B 7 holds. This implies that both / and R~j have 
bounded support. By Lemma 12.11 we have that — £ is a subordinator with a strictly positive drift 
d. Since B 7 is supported by [0, 1] and the support of / is given by [0, d^ 1 ] we have also that 
sup {t > : ¥(R > t) > 0} = d. Besides, using the expression (|17p for the moments of / and that 

E[B ? l = r 'l + f'"^) ,forn>l, 
' r(n + 7+l) 

it is also obvious that the entire moments of i? 7 are then given by (|39p . These moments characterize 
the distribution of i? 7 since it has a bounded support. 

Next, note that E[(R y /d) n+1 ] < E[(i? 7 /d) n ], for n > 1, which leads to 

0(n + l) <d 
n + 1 + 7 ~ 

and then q + Jq°°(1 — e~^ n+1 ^ x )H(dx) < ^7 for all n > 1. Letting n -4- 00, we get 11(0, 00) + (7 < ^7. 

• It remains to prove that the factorization holds when — £ is a subordinator with drift d > 
and with a Levy measure and killing rate such that 11(0, 00) + q < dj. In the proof of Theorem 11.51 
we proved the existence of a random variable, say R, independent of I such that 



RI L ^B 



(n(0,oo)+g)/d- 
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To get a similar factorization of B 7 , for 7 > n (°'^ : ')+'? =: ^, 0) note that if B 7o+ i )7 _ 7o is a Beta 
random variable independent of B 70 then 

- D 70+l,7-70- D 70 ~~ ■ D 7 - 

So, take a version of B 70+ i )7 _ 70 which is independent of R and / and set i? 7 = i?B 70+ i >7 _ 70 . We 
indeed have R y I L = w B 7 with R 1 independent of /. □ 

To finish this section we mention that for brevity further details about the random variables 
R~ will be provided somewhere else. 



5 Yaglom limits: Frechet cases and factorization of Pareto distri- 
butions 

We now turn to the Frechet cases, when the pssMp X is not monotone. 

Proof of Theorem \1.(A The assertion: / G MDAprcchot if and only if t 1— > P(J > t) is regularly 
varying at infinity with some index —7 < 0, follows from the assertion (iii) in Theorem 12.61 The 
convergence in (|15p follows from (iii) in Proposition 12.71 where the factorization of the Pareto 
distribution in terms of I is also proved. The sufficient conditions (i) and (ii) in Theorem 11.61 are 
proved respectively in [33] and [3D]. 

We are just left to prove that a necessary condition for / G MDAprechet is that Efe 7 ^ 1 ] < 1 and 
]E[ e (7+<5)€i] > 1 for all 5 > 0, where 7 > is such that —7 is the index of regular variation for the 
right tail distribution of /. If J G MDAprecheti the regular variation of the tail distribution of I 
implies that I has positive moments of all orders /3 < 7. From Lemma 3 in [30] it follows that this 
happens if and only if 

Efe^ 1 ] < 1, V/? < 7. 

Then by tacking limit as /3 f 7, and by a combination of the monotone and the dominated conver- 
gence theorems we get that 



1 > limE [e^ 1 ] = limE 

/3f7 L J /3t7 



e*l{6«)} 



+ limE 

/Sf7 



Efe^ 1 ]. 



(40) 



Also we have that E[/ 7+<5 ] = 00, for all 5 > 0, and hence Lemma 3 in [30] implies that Efe^ 7 " 1 " 15 ^ 1 ] > 1, 
and actually E[e^ 7+<5 ^ 1 ] > 1 by strict convexity, for all 5 > 0. □ 

As we did it in the other cases, we provide a necessary and sufficient condition for a factorization 
of the Pareto's distribution in terms of /, and a description of the law of J 7 in the factorization 



Proposition 5.1. For 7 > 0, there exists a unique in law random variable J-y independent of I, 
such that 



t j Law ^ 



ijff E[exp(7^i)] < 1. Furthermore, 

(i) i/Efe^ 1 ] < 1, then E[H] < 00 and 



ry f°° A X 

^ € dy) = E\lT} J ^ P *(*iedy,KT ), y>0, 
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(ii) i/E[e^] = 1, then 

P(J 7 E dy) = E 1 ^-^((J*)- 1 E dy), y > 0, 

where I* = exp(— £g)ds, where £,* is a Levy process that drifts towards +oo with distribu- 
tion P* o where P* is the probability measure defined by 

= e**P| ft , & = <7(f s ,s <t), *>0. 

The proof of this result will be given in the next section. This will be done by establishing 
a connection with quasi-stationary measures for the process of the Ornstein-Uhlenbeck type Ut = 
e~'X e f_i, which is of interest in itself. In fact we will prove that the conditions in the latter 
proposition are equivalent to the existence of a quasi-stationary measure for U. 

6 Connections with quasi- stationary distributions for Ornstein- 
Uhlenbeck type processes and proof of Proposition 15.11 

Assume that £ is a real valued Levy process drifting towards — oo, and that £ is not the negative of 
a subordinator. Recall that the index of self-similarity of the associated pssMp X is assumed to be 
1. We have seen that the existence of Yaglom limits for pssMp in the non-monotone case is closely 
related to the existence of a random variable J 7 such that 

t j Lew -p 

for some 7 > 0. In this section we would like to prove Proposition 15. 1\ that is give a necessary and 
sufficient condition for the existence of such a factor J 7 and, in the case where it exists, characterize 
its law. To that end we start by studying the problem from the perspective of Ornstein-Uhlenbeck 
type processes. 

We assume To < 00 a.s. and consider the process of the Ornstein-Uhlenbeck type associated to 
Xbj 

U t = e~ t X e t_ 1 , t>0. 
Observe that U is a process that hits in a finite time Tq . Lamperti's transformation implies that 

C 



We would like to know under which conditions there exists a quasi-stationary measure for U, viz. 
a probability measure on (0, 00), say u, and an index 9 > 0, such that 

f u(dx)E x [f(U t ),t <T C/ ] =e- et [ u(dx)f(x), t > 0, 

for any / continuous and bounded. Remark that if there exists a quasi-stationary measure associ- 
ated to U, say u, with index 8, then we have that 



e~ et 



[ v{dx)¥ x {TU >t) 
Jm* + 

/ u(dx)¥(xl > e* - 1), t>0. 
Jr* 
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It follows that if we let Jg be a random variable independent of I and with distribution u, we have 
that 

F(J e I>t) = {l + t)- e , t>0. 

That is 

So we see that there exists a random variable Jg independent of / such that IJg follows a Pareto 
distribution whenever there is a quasi-stationary measure for U. This indicates us that to tackle 
the problem of finding NASC for the existence of a factorization of the Pareto law in terms of / we 
should look for NASC for the existence of a quasi-stationary law for U. For that end we point out 
a further connection with the so-called self-similar entrance laws for a pssMp. 

We will say that a family of sigma-finite measures on (0, oo), {ijt, t > 0}, is a self-similar entrance 
law for the semigroup {P t x ,t > 0} of X if the following are satisfied 



VsP t x = m+s 



(EL-i) the identity between measures 
that is 

/ Vs (dx)E x [f(X t ),t < T ] = [ m+s (dx)f(x) 
Jr* + Jr* + 

V/ : (0, oo) — > R positive measurable, holds for any t, s > 0. 
(EL-ii) there exists an index 7 > such that 

Vsf = s'^mHsf, 

where / denotes any positive and measurable function and for c > 0, H c denotes the dilation 
operator H c f(x) = f(cx). 

(EL-iii) r/i is a probability measure. 

In that case, we say that {r/ s , s > 0}, is a 7-self-similar entrance law associated to X. We have the 
following lemma that relates the QS-measures for U and the self-similar entrance laws for X. 

Lemma 6.1. There is a bijection between the family of self-similar entrance laws associated to X 
and the quasi- stationary laws associated to U. More precisely, let v be a 7- quasi- stationary measure 
for U, and define a family of measures {i] s , s > 0} by 

n s f := s~~<vH s f. 

Then {n s ,s > 0} forms a 7- self- similar entrance law for X. Reciprocally, if {n s ,s > 0} is a ^-self- 
similar entrance law for X then v(dx) = r]i(dx) is a 7- quasi- stationary law for U. 

Proof. Let {n s , s > 0}, be a 7-self-similar entrance law associated to X. We claim that the measure 
v(dx) := 7]i(dx),x > 0, is a 7-quasi-stationary distribution for (the semigroup {P/ 7 ,t > 0} of) the 
process U. Observe that for any function / bounded and measurable we have that 

Iff = P$-iH e - t f. 
Then by the hypotheses (EL) we have that 

vFff = mP^He-tf = rj e *H e - t f = e-^ m f = e^vf. (41) 

Which proves that v is a QS measure for U. Now, let v be a 7-quasi-stationary measure for U, and 
define a family of measures {rj s , s > 0} by 

Vsf ■= s~~<vH s f. 
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We claim that rj forms an entrance law for X. Indeed, we have the following identities: 
risP* f := s-i [ v(dy)E ys [f(X t )} 

= 8-1 j u(dy)E y [f(sX t/s )] 

= s-t / u(dy) E y [/((« + t)e~ l ^ 1+t ^X eM 
it* L 

u (^y) p Hi+t/s) H (t+s)f 



a(l+t/s)_i, 



(42) 



s -'(l + t/s)-^H {t+s) f 

Vt+sf- 



□ 



Therefore, in order to characterize the quasi-stationary laws for U we need to characterize the 
self-similar entrance laws of X. This is a problem that has been studied by the second author in 
[381 139] - related to the existence of recurrent extensions of the process X, and by Vuolle-Apiala in 
[44j , In those references the parameter 7 of self-similarity for the entrance laws is restricted to be 
in (0, 1) because those are the only parameters which are relevant for the existence of recurrent 
extensions. Disregarding that restriction we can deduce from those papers that there are only two 
types of self-similar entrance laws, either 77 = (rjt,t > 0) is such that 

lim jftlr/ )} =0, a > 0, 



or 



> 0, a > 0. 



In the second case, Vuolle-Apiala |44j proved that there exists a measure /j such that rjt(dy) = 
f^°n(dx)Pj X: (x,dy), t > 0, and in fact there is a 7 > 0, and a constant < c 7 < 00 such that 
fj,(dx) = c- f x~( 1+ ' y ^dx, for x > 0. A description for 77 in the first case and £ (0, 1) has been given 

in I3SH32]. 

We have now all the elements to prove Proposition 15.11 



Proof of Proposition \5.1[ Let 7 > 0, and assume that there exists a random variable J 7 independent 
of I such that 

Law 

L/7I 'y . 

Given that a Pareto random variable of index 7 has moments of order (3 6 [0, 7) it follows that 
E[/^] < 00 for all /3 G [0,7). According to Lemma 3 in [3D] we have that the latter implies that 
E[e^] < 1, for all /3 G (0,7). Ar guing as in (jlOj) we get E[e 7?1 ] < 1, which proves the direct 
implication. In order to prove the converse implication we will first prove that the laws in (i) 
and (ii) are well defined and indeed satisfy the property of being one of the factors in the Pareto 
factorization. So, assume first that £ is such that Efe 7 ^ 1 ] < 1, that X is the pssMp associated to £ 
via Lamperti's transformation and that J 7 is a random variable independent of /, and whose law 
under P is the one described in (i). From Lemma 3 in [ID] we know that E[/ 7 ] < 00. Using that 
under Pi, To has the same law as /, it is easily verified that 

f Pi (l/u < T ) 4^- = 1. 
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It follows that the measure defined in (i) is a probability measure. From the discussion following 
the proof of Lemma 16. II we know that the family of measures (rjt,t > 0) defined by 

7 f , , du 



r, t (dy) = ^_^F u (X t ed y ,t<T l 



u 



1+7 



form a 7-self-similar entrance law for X, and therefore r/i is a QS-measure for U, and hence given 
that J 7 ~ rji, and is independent of I, we get that 

j 7 L ^_ w -p 

We next assume that 7 is such that the Cramer's condition Eje 7 ^ 1 ] = 1 is satisfied. We consider 
the Levy process whose law is P* o and P* is as defined in the statement of Proposition [5TTJ 
We denote by X* the pssMp associated to £* via Lamperti's transformation. By the optimal 
sampling theorem it follows that the absolute continuity property between P* and P. is preserved 
under Lamperti's transformation, in the sense that 

P* = X? P x , over a(X s ,s < t), for t > 0, x > 0. 

It has been proved in [41] that the measures defined by 



i4f ■= E* 



/ 



11 
1* J 1* 



for / bounded measurable, 



form an entrance law for the semigroup of (X* ,¥*). It follows from the absolute continuity between 
the semigroups of X and X* that the family of measures (r)t,t > 0) 



Vt{dx)f(x) = c 7 r 7 E* 



( r* ) (i* 



f bounded and measurable, (43) 



with c 7 a normalizing constant, form an entrance law for X, viz. 

f Vt (dx) E x [f(X s ), s<T }= [ Vt+S (dx)f(x), 

for any / bounded and measurable function. We choose c 7 = 1/E* [(J*) 7 " 1 ] , so that r\\ is a 
probability measure, this is indeed possible because if 7 > 1 then 

E *[ e -(7-i)«] =E[e ?1 ] < 1, 

and then Lemma 3 in [30] implies that E* [(/*) 7 ~ 1 ] < 00; when 7 < 1, Lemma 2 in [39] ensures 
that E* [(/*) 7_1 ] < 00. Furthermore, by construction it is plain that {f]t,t > 0} is a 7-self-similar 
entrance law. It follows from the Lemma 16.11 that 7/1 is 7-quasi-stationary law for the Ornstein- 
Uhlenbeck process associated to X, and thus, from the discussion before Lemma 16. 11 that tacking 
J 7 as a random variable independent of / and with law 771, we have U 7 L = v P 7 . It is worth 
mentioning that the latter constructed entrance law coincides with the one constructed in [39] 
under the assumption 7 € (0, 1), but the method of proof used there can not be directly extended 
to deal with the case 7 > 1. 

Finally, the uniqueness in law of J 7 follows as in the proof of Lemma 12.91 □ 



From the latter theorem and its proof we infer the following corollary. 

Corollary 6.2. For 7 > 0, there exists a 7- quasi- stationary distribution for the Ornstein-Uhlenbeck 
type process U associated to X if and only if E[exp(7^i)] < 1. In that case the 7- quasi-stationary 
distribution is the measure described in Proposition \5. 11 
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Corollary 6.3. Let the pssMp X be not monotone and U = (Ut = e~ t X e t-i,t > 0) the process of 
the Ornstein- Uhlenbeck type associated to X. We have that U admits a Yaglom limit if and only if 
I € MDApr&het! that is if and only if t ^ P(7 > t) is regularly varying at oo, say with index —7, 
7 > 0. In such a case, for x > 

¥ x (U t e- I t<T%) >^~>\ (44) 

where fij is as described in Proposition \5.1[ Necessary and sufficient conditions for I € MDAprechet 
are given in Theorem \l.b\ 



7 Examples 

In this section, some examples illustrating our results on Yaglom limits of pssMp are detailed. 
In some cases, they also lead to new (to our knowledge) factorizations of the Beta and Pareto 
distributions. 

Throughout the section X will denote a pssMp with self-similarity index 1/a > 0, and £ the 
Levy process associated to it via Lamperti's transformation. Recall that in this case X a is the 
1-pssMp associated to the Levy process a£, and hence our results can be easily translated to deal 
with the general case. 



7.1 Gumbel cases and factorizations of the exponential law 

In [22J some examples have been provided under the assumption that the underlying Levy process 
is the negative of a subordinator with infinite lifetime and a regularly varying tail Levy measure. 
Further examples can be deduced from recent papers where the law of the exponential functional of 
a subordinator and that of its factor R are determined, see for instance [5], [30], [33], [21]. Instead 
of listing those examples we will restrict ourselves to provide an example where the underlying 
subordinator has a finite lifetime. 

Example 1. Assume that under P x , X is the negative of an a-stable subordinator issued from 
x > 0, and killed at its first passage time below 0; for < a < 1. This is a 1/a-pssMp with 
infinitesimal generator given by 

L m = [ (/(* - v) - m) c -0 - £»m 



x~ a 



(jf (/(»-) -/m){£S£=-£/m) . 



for / : [0, 00) — > K. measurable, /(0) = 0, smooth enough and that vanishes at infinity. The second 
equality follows from the first by a change of variables, while the first one is obtained from the 
known expression of the infinitesimal generator of X. The underlying Levy process £ is the negative 
of a subordinator killed at an exponential time with parameter and its Laplace exponent is 
given by 

c+ r(A + i)r(i-q) 

<t> (A) = . rr-, A > 0. 

a 1 (A + (1 — a)) 

For notational convenience we chose c+ = a/r(l — a). It can be deduced from a classical result by 
Bingham |10] , Proposition 1, that under Pi the first passage time below for X has a Mittag-Leffler 
distribution with parameter a, and hence that 



C 



Law 

'a 
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where r Q follows a strictly stable law of parameter a. It has been proved by Shanbhag and Sree- 
hari [32] that the equality in law 



f\ <~v LcLW 

TVv 6 = 



holds, with e an exponential random variable independent of r a . See [5] for other extensions of this 
factorization. We have the following corollary to Theorem 11.21 

Corollary 7.1. Assume that under ¥ x , X is the negative of an a-stable subordinator issued from 
x > 0, and killed at its first passage time below 0, < a < 1. Then the tail distribution of the Levy 
measure of — £ is regularly varying at with index —a and hence Ml}) is satisfied. Therefore: 

Pi ( (c + a a )~^=^t~^X t e • I t < T ) ► P(e € •) 

7.2 Weibull cases and factorizations of the Beta (1,7) r.v. 

We just provide here an explicit example of a factorization of the Beta random variable with 
parameters (1,7), 7 > 0, further examples and properties of the factor p\ will be provided in 



Example 2. Assume that under P, — £ is a subordinator with Laplace exponent 

(p(X) = cA + = cA + / (1 - e- Xx )qpe- px dx, A > 0, (45) 

A + P Jo 

where c, q > and p > (the second equality having sense only when p > 0). Hence when p > 0, 
the Levy measure is qpe~ px , x > 0, and has total mass q, while when p = 0, — £ has no jumps and 
is killed at rate g. For a > 0, we will denote c a := ca and p a := p/a. Observe that the Laplace 
exponent of — at; is given by cp(aX), A > 0, and so it can be expressed as in the rightmost term in 
equation (|4"5|) with c and p replaced by c a and p a respectively. It is known, and easy to show via 
an identification of the moments of /, that in this case / = J °° exp(a£ s )ds has a density 

k(x) = r X Pa (l - CaX)<* l {a . e(0) ( Ca) -i) } , ICR, 



r (ft a + 1) r 



that is / has the same law as (c a ) 1 B pa+1 _q_, see e.g. [13]. According to Proposition 14.21 for any 
7 > there exists a r.v. i? 7 such that 



and the moments of R-y are given by 



(^Kl = W-"IIrS=II 



L aw 



n 1 / .\ n 1 1 o 



2=1 1=1 



r(n + i) r(i + 7 ) r(n + Pa + ^ + i) r(i + Pa ) 
r(i) r(n + 7 + i) r( Pa + ^- + i) r(n + Pa + i) 

(l)n (4 



(46) 



(l+ 7 )n(^ + d)n : 



where (/3) s := pTgy) denotes the Pochhammer's symbol, and £ := p a + 1 + ^— , (f = — We deal 
first with the case p Q > 0. These parameters obviously satisfy that £, 7 + d > 0, 1 + 7, 1 + p a > 0, 
and that min{l,/3 a + 1 + -?-} < min{l + 7, 1 + p a }. This simple remark allows to ensure that 
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corresponds to the moments of a BetaProd random 
r-), that we will denote by B 1 „ . + _q_ , 1 q_. BetaProd 



the expression of the rightmost term in 
variable with parameters (l,7,p Q , + ^- + l, 

random variables were introduced by Dufresne in [19] and he proved that these random variables 
are determined by their entire moments, calculated its Mellin's transform and calculated explicitly 
its density in terms of the Gauss hypergeometric function, which we do not reproduce here. It 

As a consequence of Proposition 14.21 we get the identity 



follows that R, 



Law 



1,7, Pa + 



-l,^- B l, 7 ,p Q + ^-+l - JL 



Law 



B^y. 



correspond to those of a Beta random 
c a B 1+ _2_ 7 __2_ . Finally, in the case where 



i 



In the case where p = 0, and 7 > — , the moments in 

Cot 

variable with parameters (1 + 7 — and hence R-y 

Cot Cot ' 

p = 0, and 7 = -2-, the moments in (|46p are all equal to 1, and hence i? 7 = c 

We have all the elements to state the following corollary which is an easy consequence of 
Theorem 11.51 

Corollary 7.2. Let X be the a-pssMp associated to £, whose Laplace exponent is given by 4> as in 
equation We have that X has the following Yaglom limit 



( 



lim 

t-»l/c 



X t 



--t 



l/a 



t < Tn 



3 



l/a 



7.3 Prechet cases and factorizations of Pareto distributions 

Here we include just a few examples of pssMp that have found applications in other areas and/or for 
which it is possible to determine explicitly the law of the first hitting time of zero and of its Yaglom 
limit law, and hence that give place to an explicit factorization of Pareto r.v. Other examples will 
be given in [23] . Further examples can be extracted from recent literature where an important 
effort has been made to obtain explicit distributional properties of exponential functionals of Levy 
processes, see for instance [8], [25], [32], [33] and the reference therein. 



Example 3. Let £ = (crBf — bt, t > 0) with B a Brownian motion, a 7^ 0, b > 0, and for a > 0, 
X the 1/a-pssMp associated to £ via Lamperti's transformation, i.e. a R+ diffusion process with 
infinitesimal generator: 



l f(x) = y x f ( x ) + ( y 



b ) x l - a f{x) 



x > 0, 



(47) 



for / : [0, 00) — > K smooth enough, /(0) = 0. It is well known that if < 6 < 1, a = 1 and a = 2, X 
is a 2(1 — 6)-dimensional Bessel process, and in particular when b = 1/2, X is a Brownian motion 
killed at 0. The associated Ornstein-Uhlenbeck type process Ut '■= e~ l l a X e t_x has an infinitesimal 
generator 

.l-a 



L u m = ^-y"(x) + ((^ 



- ) /'(*), 

a 



x > 0. 



(48) 



We have that the conditions in Theorem 11.61 (i) are satisfied by a£ with 7 = 2b/(aa 2 ), and hence 

= (a(o~Bt + bt),t > 0) (see Proposition 15.11 for the definition of notation £*, J*). Dufresne 
established that the exponential functionals 

poo poo 

1=1 exp(a(aB t - bt))dt, I* = exp(-a(aB t + bt))dt, 
Jo Jo 
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both have the same distribution as -, — vr-^ , with 'yoh/tn^ that follows a Gamma distribution 

with parameters (2&/(a<7 2 ), 1). The random variable ^/(ao- 2 ) defined in Proposition 15.11 is then 
such that 



(ao? 



and we have the identity in law 

1 Law ---j 

e = ^2b/(aa 2 )- 

726/( QCT 2) 

Corollary 7.3. Let X be the diffusion process killed at its first hitting time o/O, whose infinitesimal 
generator is given by (f^7| ). We have that for any x > 0, 

\ E dy | t < T ) fi(dy) := "r^^Vo}*, c = (a<r) 2 /2. 



TTie law // defined above is a quasi- stationary law for the Ornstein- Uhlenbeck type process U with 
infinitesimal generator f^g| ) and for any x > 0, 

f rj E I + ^ ^ 



i— »oo 



Observe that this result applies to the totality of self-similar diffusions that hit in a finite 
time because the only Levy processes with continuous paths are Brownian motion with drift. The 
result in the second part of the above Corollary is closely related to a seminal result by Mandl [29] 
for the classical Ornstein-Uhlenbeck process. 

Example 4 (Stable Continuous state Branching processes). Let X be an a-stable continuous state 
branching process, with a £ (1,2), that is a strong Markov process with values in [0, oo) having 
the branching property and hence with Laplace transform 



E, 



e -xx t 



exp(— xut(X)), A > 0, 



and ut(X) determined by the equation 

1 

-du = t, t > 0, 



where c+ > 0, is a constant. Kyprianou and Pardo [26] and Patie [31] showed that X is a positive 
self-similar Markov process with self-similarity index l/(a— 1), and that its underlying Levy process 
£ has no negative jumps and a Laplace transform given by 



E 



e~^ 



T(A + a) , 



where E[— £i] = m = c+r(a)T(— a) > and by convention l/r(0) = l/r(— 1) = 0. Observe that in 
this case the Levy process £* of Proposition 15. II is also spectrally positive and its Laplace transform 
is given by 



E [exp(-Ag)] = E [exp(-(A - 1)&)] = exp ( tm^ X + " ^ 



r(A- i)r(a) 



. r(A + a — 1) . 
exp ( tm(X - 1) r(A)r(a) ) , A>0, 



(49) 



In [26] Lemma 1 and in [HI] it has been proved that the first hitting time of for X, Tq L ^ 
I = J °°exp((a — l)£ s )ds is such that (c+(q — 1))/ follows a Frechet distribution with parameter 



35 



l/(a — 1). They furthermore proved that X admits a Yaglom limit, more precisely they obtained 
that for x > 0, 



lim E a 

i— >oo 



( 1 + A -(«-i))V(«-i)' 



A > 0. 



In [23] it has been proved that the rightmost term in the above display is the Laplace transform of 
a r.v. £( a _i), whose tail distribution is given by 

P (S (Q _ 1} > S j - ^ - s > 0. 

Theorem 11.61 gives another way to obtain this Yaglom limit, observing that X a ~ l is the 1-pssMp 
associated to the Levy process (a — 1)£, and that from the expression of the Laplace exponent of 
£, the Cramer's condition and the integrability condition in Theorem II, 61 - (i) are satisfied by £ with 
index 7 = 1, and hence by (a — 1)£ with index 7' = l/(a — 1). With this approach, the weak limit 

of the law of X?' 1 /t \ t < To is that denoted by p,j 1/c< * x) . An alternative characterization of this 



limiting distribution is given in 

We can deduce from this and the discussion above the following factorization of the Pareto 

random variable. Since pf 1 ^"^ 1 ^ is the law of c+(a — l)S?^j L 1 x, we have that when Wi/(q,_i) 
follows a Frechet distribution with parameter l/(a — 1) and it is independent of X^.^, then 

^i/(-i) s r-i) L = Wp ^i- 

Example 5 (Stable process killed at (— 00, 0)). Let Y be an a-stable Levy process, < a < 2, 
with positivity parameter p := P(Yi > 0), hence the Levy measure of Y is 

dx dx 

with c+, c_ > 0, c+ + c_ > 0. Let X be the process obtained by killing Y at its first hitting time 
of ( — oo,0), say Tr—oofi)- X inherits the scaling property and strong Markov property from Y and 
hence is a pssMp with self-similarity index 1/a. Chaumont and Caballero [11] established that the 
Levy process £ associated to X via Lamperti's transformation has Levy measure 

and killing rate c_/a. It can be deduced from the results by Caballero et al. [12], who introduced 
the so called Lamperti-stable class and of which this process is an element, that the characteristic 
exponent of £ takes the form 

E[ ^ 11 . exp ^ r( _ Q) E(^ + c . r( _ a) _Q^), AeK . (50) 

A similar expression of the exponent can be read in the paper [23] , where a particular choice of the 
constants c + and c_ is made. Our results can be applied to the 1-pssMp X a . Indeed, in [11] it has 
been proved that the hypotheses (i) in Theorem 11.61 are satisfied by £ with 7 = a(l — p), and hence 
by a£ with 7 = (1 — p). We obtain from Theorem 11.61 that X a admits a Yaglom limit, and from 
Proposition 15.11 that the limit law 771 is given by 

mif) = E[(I*)-p] ^/(VnUTa / > 0, measurable. 
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From [3S] Example 3, we know that the measure r/i defined by 

W) ■= E[(J L p] nf(l/(I*) 1/a )(lT P ], / > 0, measurable, (51) 

is the law of the excursion process, associated to the excursions of Y from its past infimum, at time 
one conditioned to have a lifetime larger than 1. Said otherwise, the law r\\ is the law of the stable 
meander of length one at time one, see e.g. [H] for further details about the meander process. Since 
the processes Y and X coincide before To, the above discussion establishes the following corollary. 

Corollary 7.4. Let Y be an a-stable Levy process, with < a < 2 and positivity parameter 
p = P(Yi > 0), T^oo Q) the first passage time below for Y, and Z\ be the a-stable meander process 
of length one at time one. We have the following convergence, for x > 

^ G • | t < T^o)) P(Zi G ■)■ (52) 

Furthermore, the following factorization of the Pareto distribution holds 

T^ 0) Z? L = Pi-,, (53) 

where T}^ q ^ denotes the first passage time below for Y issued from 1 and the factors on the 
left-hand side of the equality are assumed independent. 

The result in (|52p is a particular case of Lemma 15 in [T7]. In the case where Y has no negative 
or no positive jumps it is possible to obtain further information about the law of the a-stable 
meander of length one at time one. For brevity we omit the details and refer to |23j. 

Example 6. Let £ be a Levy process with drift +1 and no-positive jumps, and Levy measure 

n(-oo,-x) = be- x{b - s \ x > 0, 
with < 5 < b. The Laplace exponent of £ is 

In (E[exp(A6)]) = ^(A) = A > 0. 

The conditions in Theorem ll.61 (i) are satisfied with y = 5. The Levy process ^* has Laplace 
exponent 

ln(E[exp(A£)]) = ^*(A) = A > 0. 

According to |43j the exponential functional / has the law 

!-.= r e ^d S L = awi " B ^ +i 



/0 Rs,b-5+l 

While using the results in [7] it can be easily verified that 

r := fV«d s L 



'o &6,b-5 
and hence the random variable J$ from (ii) in Proposition 15.11 is such that 

T Law -j-. 

J& = n\,b-6- 
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Corollary 7.5. Let X be the 1-pssMp associated via Lamperti's transformation to the spectrally 
negative Levy process £ with Laplace exponent tp. For x > 0, we have 



x \ t 



^edy\t<T ) ► P (Bi )6 _i € dy) 



t—>oo 



Furthermore, the following factorizations hold 



1 — &5,b-6+ 1 r, Law 
— — -_ &l,b-6 = *5- 



5,6-5+1 
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